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A NUMERICAL ANALYSIS OF THE USE OF PERFORATED WALLS 
TO CONTROL SHOCK LOCATION AND MOVEMENT IN AN 
INTERNAL- COMPRESSION SUPERSONIC INLET 

By Clarence W. Matthews 
Langley Research Center 

SUMMARY 

This paper presents an analysis of the shock motions in an internal- compression 
supersonic inlet with the walls in the supersonic section so designed that they are perfo- 
rated to subsonic flow and solid to supersonic flow. The shock motion is computed with 
a computer program utilizing one-dimensional unsteady-flow characteristic theory. The 
perforated walls were studied to determine their ability to control the shock motion, and 
especially to contain the shock within the inlet for the application of various flow tran- 
sients to the engine face. 

It was found that the perforated walls controlled the shock motion, thereby con- 
taining the shock within the inlet for much larger transient magnitudes than was possible 
with solid walls. The response of the shock motion to the engine-face transients indi- 
cated that the response characteristics for short impulse-type transients were sufficiently 
linear to allow the use of a transfer function over a limited range. However, period dif- 
ferences and the development of a limit cycle in the shock motion when the shock position 
is ahead of the throat seriously restrict or even prevent the use of transfer functions to 
express the shock motion for long periods or throttle-chop transients. 

INTRODUCTION 

Supersonic inlets which incorporate large amounts of internal contraction operate 
most efficiently if the normal shock is maintained at the throat. However, if the shock 
moves upstream of the throat, the flow in the entire inlet breaks down. Hence, control 
of the normal- shock system to prevent breakdown of the flow in the inlet is imperative if 
the inlet is to operate efficiently. Various methods of controlling this shock motion have 
been studied experimentally, such as the use of pressure sensors or probes or the use of 
subsonic Mach numbers to vary either the bypass flow or the local area distribution of 
the inlet (refs. 1 to 5). Another method of control which avoids elaborate mechanisms is 
the use of aerodynamically variable or ''educated" holes in the walls of the converging 



supersonic portion of the inlet which bleed much of the flow in the subsonic portion 
behind the shock but bleed little of the supersonic flow ahead of the shock (refs. 6 and 7). 
An advantage of this method is that it acts as fast as the shock motion itself and so has 
the potentiality of preventing possible unstarts of the inlet in case severe engine tran- 
sients occur. Also, after such a transient ceases and operation returns to normal, the 
"perforated” wall has the potentiality of quickly returning the inlet flow to its normal 
operating condition. 

This type of inlet was studied in some detail in the experiments of reference 7, 
which gives pressure recovery results, ability to restart, and in several cases the ability 
to hold the shock stable in the contracting portion of the diffuser. It does not, however, 
present much information on the actual shock motion due to various engine transients or 
on the increase in the amplitude of such transients that the wall is capable of controlling. 

It is the purpose of this paper to present a study of the use of perforated walls in 
an internal-compression inlet to control the motion of the normal shock due to the appli- 
cation of various transient flow-property changes at the engine face. The analysis con- 
siders various degrees of wall perforation up to 40 percent open, as well as transient 
changes in velocity, pressure, and mass flow of as much as 50 to 60 percent at the engine 
face. Some discussion of the computer program used to compile the unsteady flow field 
is included. 


SYMBOLS 

A local cross-sectional area of inlet, meters 

A c h restricted area at engine face such that flow is choked 

a nondimensional sonic speed (nondimensionalized with respect to sonic speed 

at initial free-stream conditions) 

c speed of sound, meters/sec 

d diameter of area equivalent to cross section of inlet, meters 

F body force, newtons 

DS 

f entropy function, — = f(a,u,S, £,t) 

Dt 

h height of entrance to inlet, meters 


f 


/ 


2 



k coefficient having the sign of the transient disturbance and proportional 

to its magnitude 

M Mach number 

m e mass flow rate at engine face nondimensionalized with respect to mass 

flow rate at entrance of inlet 

m w mass flow rate per unit length of perforated upper wall of inlet, 

kilograms / meter- second 

P = era + u 

p static pressure, newtons/meter^ 

p e static pressure at engine face nondimensionalized with respect to free- 

stream total pressure 


Q = oa - u 


u 


nondimensional value of entropy, 


Spe ci fic entropy 


(y - 1) (Specific heat at constant pressure) 


nondimensional local velocity (nondimensionalized with respect to sonic 
speed at initial free-stream conditions) 


Uc 


value of u at engine face 


W 


nondimensional shock velocity (nondimensionalized with respect to sonic 
speed at initial free-stream conditions) 


w 


width of inlet, meters 


length coordinate along inlet axis, meters 


r 

6t 

6_ 

67- 


ratio of specific heats 

line derivative along P-line (a line of slope u + a), + (u + a)-|^ 


'H 


3 8 

line derivative along Q-line (a line of slope u - a), -^ + (u - a)-^ 


3 



nondimensional substantial derivative of entropy, that is, derivative along a 

streamline, — + u 

’ 9t a| 

density of fluid, kilograms/meter^ 

nondimensional time; unit value of r equals h/c 

SI' nondimensional value of mass removed per unit length in unit time, 

hc oo m w/v A Poo 

i// local proportion of open area to total area of walls 

u> angular frequency of sine or cosine transients 

Subscripts: 

f in front of shock 

h value of flow property in wall perforation 

r at rear of shock 

s at shock 

°° free- stream condition 


DS 
D T 



T 


THEORY 

Perforated Walls for Control of Shock Motion 

Control of the shock motion in an internal- compression supersonic inlet is neces- 
sary because serious losses in mass flow, pressure recovery, and drag effects occur if 
the shock is expelled from the inlet. It is well known that such shock expulsion does 
occur whenever the shock is forced ahead of the throat of the inlet by either engine-face 
or inlet- entrance transients in the flow. This expulsion occurs primarily because the 
shock position is unstable in a region wherein the derivative of the area with respect to 
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distance along the inlet is negative. A possible means of preventing this expulsion is to 
reduce the pressure behind the shock, thereby drawing the shock back to its stable region 
downstream of the throat. This result may be accomplished by bleeding the air from the 
inlet through perforations in the wall. 

References 6 and 7 present results of studies of the behavior of inlets with perfo- 
rated walls. These studies, which were mostly experimental, indicate that the shock will 
remain within the inlet whenever a transient such as a throttle chop is applied to the inlet. 
However, because of experimental limitations, no detailed study of the actual shock 
motion caused by various transients has been made. Such an analysis would be difficult 
to make experimentally, but can be accomplished by computation with a one-dimensional 
time-variant characteristic system. Although several of the studies used an inlet with 
"educated" holes which bled much of the subsonic flow but little of the supersonic exit 
flow, the present analysis considers an inlet which has spring-loaded flap doors over the 
surface of the inlet. (See figs. 1(a) and 1(b).) These doors are considered to have suf- 
ficient spring loading to close them when the low-pressure supersonic flow covers the 
doors and to open them when the high-pressure subsonic flow covers the doors. Although 
the physical realization of such a door might present many problems, such as unstable 
flapping or removal of the fluid which escapes through the walls, the condition is simply 
expressed analytically by setting \f/, the proportion of open area in the wall, equal to 
zero if the flow is supersonic and equal to some given value if the flow is subsonic. It 
is assumed that the pressure exterior to the inlet is small enough to allow any flow 
through an opening to be sonic, and that the holes are completely filled with flow. 

Computing Program 

The program used to compute the shock motion in an inlet due to various distur- 
bances involves an integration over a £,t field, where £ is the nondimensional dis- 
tance along the axis of the inlet and t is the nondimensional time, of the one- 
dimensional unsteady-flow characteristic equations presented in reference 8. These 
equations are expressed in terms of the sonic speed a, the flow velocity u, and the 
entropy S. The values u and a are nondimensionalized with respect to the speed of 
sound at free-stream conditions, that is, a = 1 at initial free-stream conditions. The 
unit distance for £ is considered to be the height of the two-dimensional inlet at the 
inlet entrance. The corresponding unit value of r is the time for a free-stream sound 
wave to traverse the unit distance for £. The characteristic equations actually express 
the changes in P and Q, functions of u and a, along the two characteristic directions 

^P- = u + a and = u - a, and the change in S along the streamline P- = u. 
dr dr D dr 
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Flow Equations 


Equations for chara cteristic phenomena.- The unsteady one -dimensional flow 
equations may be integrated by a characteristic method similar to that used for two- 
dimensional or three-dimensional rotational symmetric flow. The development of the 
characteristic method and the equations used are given in reference 8. These character- 
istic equations are: 


M 

5t 


au 8A 
A 


• dr 


a 

A 


a M + (r - i )a e| +F _ 


’S'a 


(y-3)/(y-l) e yS 


(1) 


and 


5-Q _ au 9A 
6 t ‘ A 


M + a — + (y - l)a 
A 8t 5t Ky ' D t 


^a 


(y-3)/(y-l) e yS 


( 2 ) 


where P = era + u, Q = era - u, and F is a body force acting on the fluid. Equation (1) 
applies along the P-line, whose slope is u + a, and equation (2) applies along the Q-line, 
whose slope is u - a. These lines are shown in figure 2, which presents a diagram of 
the basic characteristic net. The entropy equation 

|| = f(a,u,S,|,T) 

applies along the streamline u and is given by the problem conditions. 

The term 'I' is defined as the nondimensional value of mass removed per unit 
length in unit time. This term is related to the geometry and flow characteristics of the 
inlet in the following fashion. Consider the evaluation of ^ for the two-dimensional 
inlet shown in figure 1(c). (A similar consideration could be applied to a circular inlet.) 
The mass flow rate through unit length of the perforated upper wall of this inlet may be 
expressed as 

m w = ^Ph c h 

where 


xf/ ratio of open area in wall to total area 

w width of inlet (see fig. 1(c)) 

density of fluid in wall opening 
c, local sound speed in wall opening 
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In the equation for m w , is equal to y® where c*, is the free-stream 

sonic speed; hence 

m w = ^P h \ C oo 

It is shown on page 15 of reference 8, that 

p h - < 5 » a h < ’ e ‘ ySh 


where is the density for the free-stream condition wherein a^ = 1. Thus 

i ~ ySu 

m w = i/AvPo 0 Co O a h CT+1 e h 

The term SI> is given on page 6 of reference 8. In the notation of the present 
paper, 

.p _ hC oo m W 

rAp^ 

where p^ is free-stream pressure. The term a^ is evaluated by considering that 
the flow through the holes occurs at its local sonic velocity. This value is obtained by 
first writing Bernoulli's equation in its nondimensional form, which is 

ga 2 + u 2 s Constant 

If it is considered that the flow expands so that a^ = a = u, then 

ga 2 + u 2 = ga^ 2 + a^ 2 = (a + ^a^ 2 


or 



ga 2 + u 2 
g + 1 


Substituting for m w and a^ results in 


t hc oo / ( ga 2 + u 2 \ 

^ = i °° i/zwp^c 1 ~ u 


yAp c 


(g+l)/2 


g + 1 / 


-ys h 


Note that 


^ = c 2 

Poo 


"oo 


wh = A„ 
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III 



so that 


* = \p 


, / 2 2\ (CT+1)/2 * 
Aoofaa^ + u^^ „-ySh 

A \ d+1 / e 


With the substitution of y = g * — , 

a ( r -3)/(y-l) = ^ CT " 1 

Hence the term ^a^ ^e^® of equations (1) and (2) may be expressed in terms 

of geometry and flow properties as 

Aqq /aa 2 + u 2\ ( g+1 )/ 2 i 
W A \ o+l ) a a-l 


The treatment of tunnels with other shapes or with different perforation distributions 
may change the area factor of this equation. For instance, if the inlet is assumed to be 
circular with perforations over the entire area and the entrance diameter is considered 
the unit length, then Aoo/a in the above equation must be replaced with 4 ^Aoo/A. 

DS 

The term (y - l)a — — ± F is shown on page 49 of reference 8 to contribute the 
Dr 

following to the value of A + P and A_Q, respectively: 


and 


where 


F[l - (y- 1)|a + t+ aA + S 


-F 


1 + (y - 1)^ 


A t + a A S 


F = -^u 2 

2d |u| 


Here d 



is the diameter of an area equivalent to the cross-sectional area of the inlet, or 
and X is a solution of the following equation: 

1 


ri/2 


= 2 log 10 ^ A 1/2 - 0.8 


where v is the kinematic viscosity. The terms A + S and A_S can be evaluated once 
the entropy is determined. 


The determination of the entropy proceeds in this manner. The value of the entropy 
is already known at the previously calculated points B and C of figure 2. It is desired to 
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calculate the entropy at the unknown point D. The rate of change of entropy along the 
streamline of slope u from the point T located on line BC to the point D is given by 

PS _ uF 
Dr ‘ a 2 

This value may be calculated at T by interpolations of the values at B and C and at the 
point D after a and u are estimated or calculated at that point. The value of S at 
D is then determined from the relation 


Sd = St 


uF' 


+ 



t T 


where the subscripts T and D indicate the values of the respective variables at 
points T and D. Once Sj) is known, then 

A + S = S D - S B 

and 


A_S = S D - S c 


where the subscripts B and C indicate the values of the variables at points B and C. 
After these modifications have been completed, equations (1) and (2) can be written in the 
delta form: 


A . P = - 


and 


<l( u M + M Ah 

9| a T ) h 




aa- 


2 + u 2\( ct+1 )/ 2 

<7+1 ) 


^I A + T + 


'(l-kL 


Du 


A , T + SL A.S 


( 3 ) 


A Q = -i.fu— +-^A r - ^ Aco (°* 2 + , U 2 . N ) 

A\ a§ + d T ) - A \ o+l ) 


(a+i)/2 


.<7-1 


,. T . F ( 1+ l2^!ls) 


^ ^ y A^_T + a. A_S 

( 4 ) 


The application of these equations to an actual flow field is seen upon consideration 
of figure 2. Assume that the flow properties a, u, and S are known at both of the 
points B and C. It is desired to determine the flow at point D from the known data. The 
point D may be located by the intersection of the lines of slope, u + a extending from 
point B, and u - a extending from point C. Then the equation for 6 + P/8t is integrated 
along the u + a line and the equation for 6_Q/6 t is similarly integrated along the 
u - a line, to give the new values of P and Q at the point D. A third point T is 
interpolated between B and C so that a line of slope u extending from T will intersect 
D. The entropy equation is then integrated along that line to obtain the final values of a, 
u, and S. 


i 



The integrations are actually computed by writing the differential equations in a 
delta form. Then 


P D = P B + A + P = P B + (M) b(T d - t b ) 
Qp = Q c + A_Q = Q c + (%t) C ( T D ' T c) > 

and 

S D = S T + AS = S T + - r T J 


( 5 ) 


where the subscripts refer to points in figure 2. Using the characteristic equations in 
this form permits the evaluation of P B and Q B . Once these are known, a and u 
may be determined from the simultaneous equations 

aa D + U D = P D 


CTa D " U D = Qd 
or 

„ _ P D + % 

_ p b - Qd 

U D 2 

The values of £, r, a, u, and S for point D can be determined by substituting 
the corresponding values at points B and C into equations (3) and (4) and performing the 
indicated computations. However, greater accuracy can be attained if the slope of the 
characteristic lines is made equal to the average of the slopes at each of the end points, 
so that the line slope is representative of the slopes at both B and D rather than at B as 
in the primary computation. Such a consideration of both slopes was met by computing 
each term of equations (3) and (4), using the appropriate values of A, a, u, S, and so 
forth, for each of the points B, C, and D. The appropriate values of the terms were then 
averaged to give the final values of the terms used in equations (3) and (4). The use of 
these values requires that A B , u B , S B , and be recomputed several times to 

insure stable values and thus increase the accuracy of the computation. It was found that 
about five or six such iterations would determine the value of u^ to within 0.00001. 

Equations for shock phenomena.- Normal shocks occur in the inlet as a result of 
the intersection of either two P-lines or two Q- lines. The primary normal shock which 
is necessary to satisfy the initial conditions of the inlet flow is the result of the intersec- 
tion of the supersonic and subsonic Q-lines. Other shocks may exist on either side of 
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the primary shock and are usually much weaker unless a sudden shock-type change is 
made at either the entrance or the exit of the inlet. This study is interested primarily in 
the control of the motion of the primary shock, as such motion can have serious effects 
on the flow properties of the inlet, especially if unstarting of the inlet is involved. 

The motions of either the primary shock or the subsequently formed secondary 
shocks are computed by using the equations and, to a varied extent, the methods given in 
reference 8. These equations are applied to the field shown in figure 3(a), which repre- 
sents a typical Q- shock characteristic diagram. A magnified portion of figure 3(a) is 
shown in figure 3(b). 

The discussion in reference 8 shows that P- and Q- shocks have respectively differ- 
ent values of P or Q on each side of the shock. This is expected because the two P- 
or Q- lines which intersect on the shock surface will not have the same value of P or Q. 
The shock properties, that is, Mach number, velocity, and strength, are primarily depen- 
dent on the value of A + P or A_Q which occurs across the shock. Hence it is neces- 
sary to compute the value of A + P or A_Q across the shock before its properties can 
be determined. The computation of a Q-shock proceeds somewhat in the following man- 
ner. In figure 3(b), the values of a, u, and S in the flow field at the characteristic 
points A and C are known, as well as the values at the points Bf and B r (in front of and 
at the rear of the shock). It is now desired to compute a, u, S, £, and r at the 
points Zf and Z r (the front and rear sides of the shock point Z). 

The shock point is approximately located by computing the location of the inter- 
section of the P-line from A and the shock from B, using the shock velocity W for the 
slope of the shock. Then the value of P at Zf may be computed from equation (3), and 
by using the interpolated point E whose Q-line intersects the point Zf, it is possible to 
compute the value of Q at Zf. These input values Pf and Qf are then used to deter- 
mine af and Uf. The interpolated point G may be used similarly to compute Q r , but 
P r cannot be computed. Hence it is now necessary to resort to the shock equations 
relating the shock Mach number and A_Q or A + P. The equation for AQ across a 
Q-shock may be expressed as 

AQ = Qp — Qf = cxa r — Up — CTa f ■+• Uf 

It is shown in reference 8 that 


u f = W f + a f M s>f 
u r = Wp + a r M Sj p 

where M s f is the Mach number of the flow into the shock relative to the shock and 
M s ,r is the Mach number of the flow away from the shock relative to the shock. Then 
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( 6 ) 


r + M s,f 


The substitution of a for 


into equations (95) and (96) of reference 9 results in 


7 ~ 1 

ar r^ 2 > M I,--1 1/2 ( M I,-) 1/2 
(cr + l)M g f 


af 


and 


M 


s,r 


M lf + CT 


(a + 2 )m| f 


1/2 


respectively. The substitution of these relations into equation (6) leads to 

AQ J [(a + 2)M| f - l] 1/2 (M|,f + °) 1/2 + M| >f - (o + 1)M m - l | 

(o + l)M g f 


( 7 ) 


The value of AP/af across a P-shock can be found in like manner. 

Equation (7) may be solved for M g f by using an iteration technique such as the 
Newton- Raphson. Once M„ f is known, a r may be computed and then u r may be 
determined from the shock velocity equation; that is, for a Q-shock 

u r = W + a r M g ^ r 

where M„ „ is the shock Mach number behind the shock which corresponds to M,, f in 
front of the shock and is determined from the relation 


M 


s,r 


M2+a 


11/2 


(a + 2)Mg f 


Thus, the complete properties of the shock may be determined. As with the char- 
acteristics, the various velocities may be averaged and a new shock value may be deter- 
mined. This process may be iterated to the desired degree of accuracy. The P- shocks 
may be computed in a similar fashion, that is, by first determining Pf, Q f , and P r 
from appropriate known points. Then the shock properties may be determined from cor- 
responding shock equations. However, the shock velocity relations for the P-shock are 


W = Uj + afM g j = u r + a r M s r 
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A problem that may occur, especially in computation of Q-shocks, is the problem 
of multiple intersections of the shock and the Q characteristic lines. An example is 
shown in figure 4 . In this figure the shock point Z would normally be computed from the 
known points A2, B2, and Sj, and Q r would be obtained from an interpolated point 
occurring between Si and C\. However, it is seen that the Q-lines from B3, B4, Bg, Cj, 
C2, and C3 and the shock from S2 all intersect with the shock from Si before that shock 
intersects with the characteristic line P2. Although it is possible to search out, test for, 
and calculate each separate intersection, the process would greatly complicate the pro- 
gram. This complexity may be avoided by obtaining the Qf value from interpolation 
between the points A2 and S^ The value of Q r needed to complete the shock point 
computation is obtained by checking each pair of points (S^f,Ci), (Ci,C2), (^2^3), 
(C3,S2 ? f), and (S2 ?r ,C4^ until two points are found for which a Q-line originating at Z r 
will intersect the line between the two points. Once this condition is discovered, Q r 
may be calculated and the shock values at Z may be determined. Although this method 
is not quite as exact as the method of computing each intersection, the errors involved 
should be small and involve only slight changes in shock location. Only small errors are 
expected because the values of Pf, Qf, and Q r are determined from practically the 
same characteristic lines regardless of the shock- characteristic intersections that occur 
in between. 

Input d ata for program .- The computing program allows a number of input varia- 
tions, such as time variation of velocity, pressure, or mass flow at the engine face; time 
variation of the engine-face area with choked flow assumed at that point; time variation 
of sonic speed or flow velocity at the inlet entrance; or a possible time variation of the 
area distribution of the inlet. Several time-variation functions are included. These are 
steady state, sinusoidal variation, stall, exponential ramp, linear ramp, flat-top impulse, 
and a triangular- shaped impulse. Any of these are available through a data card. The 
area derivative and the inlet wall perforation distributions are entered as tabular data. 

A cubic interpolation routine is used to interpolate the area distributions, whereas a 
linear routine is used to interpolate the perforation distribution. The linear distribution 
is used for the perforation table because the values of this distribution may change too 
abruptly for cubic interpolation. A given value of y, the ratio of specific heats, may be 
entered into the program. No provision is made for the variation of its value as this 
program is more concerned with the effects of external influences which may be used to 
control the inlet flow properties. 

Only one inlet was considered for analysis. This inlet was an internal-compression 
inlet with an area ratio reduction from 1.0 at the entrance to 0.30 at the throat. It was 
derived, with modification for friction effects, from a two-dimensional characteristic com- 
putation of a supersonic diffuser which would reduce the Mach number of the flow from 3.0 
to 1 . 3 . Several small linear area adjustments were made to the area distribution until the 
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flow could be computed without choked phenomena occurring. A linear subsonic diffuser 
with a slope of 0.08497 was added to the supersonic portion. The area distribution of this 
inlet is presented in figure 5. 

It was found necessary to fix the dimensions of the inlet if viscous effects were to 
be computed. The entrance height was assumed to be 1 meter, and the operating altitude 
was assumed to be 20 000 meters. 

Several perforation distributions were considered. A wall with no perforation was 
used for comparison. Walls with 20 percent and 40 percent perforation from the entrance 
to station 4.4 were investigated first. Analysis of steady- state flow computations with 
the shock located at various points along the inlet showed that these distributions would 
not always permit the shock to remain stationary at any point in the inlet. Hence a third 
distribution (fig. 5) was determined by trying and modifying perforation distributions until 
a distribution was attained that would allow the shock to be set at any point in the inlet. A 
shock located ahead of the inlet was treated in the following fashion. The inlet was con- 
sidered to be extended forward into the flow sufficiently to be ahead of the shock. The 
perforation of this portion of the inlet was considered to be zero if the flow was super- 
sonic and equal to 1 or fully open if the flow was subsonic. Although more sophisticated 
means might be used to represent the escaped flow ahead of the inlet, this very rough 
method has the advantage that no additional programing is necessary. Also, once the 
shock is expelled, its exact position is not important. Although some errors will occur 
in shock position and possibly in time to restart, they will not seriously affect the purpose 
of the analysis. 


RESULTS AND DISCUSSION 

Because of the nature of the characteristic programs involved, it was possible to 
select only one particular condition at the engine face and allow it to vary with time. 
Variations were applied to flow velocity, static pressure, and mass-flow rate. At the 
inlet entrance either sonic velocity or flow velocity or both could be chosen for transient 
variation. All of the studies involved in this analysis are concerned with the various 
engine-face transients. Any time variation desired could be included in the computations 
by modifying the transient subprogram. Those chosen were steady state, a sinusoidal 
variation, a stall transient, a ramp throttle change, and several impulse transients, that 
is, flat-top, triangular, and half-sine-wave. These transients cover most of those of 
interest in basic inlet control studies. 

The results of a number of computations with simulated impulse transients are 
shown in figures 6 to 10, which present shock position as a function of time for a single 
half-sine-wave pulse with a duration of 5 time units. The results of the computations for 
several stall transients and ramp transients are presented in figures 11 to 17. 
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Shock Motion Due to Half-Sine-Wave Transient 


The transient chosen for this portion of the investigation was forced on the inlet by 
allowing the engine-face condition to vary with the following equations: 


C e = C ss (l + k sin cor) 
Ce = ^ss 
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where C e is the imposed engine-face condition, C ss is the steady-state engine-face 
condition, and k has the magnitude and sign of the transient. Examination of the equa- 
tion shows that the transient has the form of a half- sine wave, and then conditions return 
to normal after the transient has passed. The engine-face values of either velocity, pres- 
sure, mass flow rate, or choking area were used to furnish the value of C ss . 

This transient function was used because it simulates the application of an impulse 
to the inlet, especially if the duration of the transient is fairly short compared with the 
duration of the response. The response to this transient therefore simulates the transfer 
function which exists between the engine-face flow property and the shock motion. 


Examination of this response can provide information on properties of the transfer 
function such as its oscillatory nature, rate of damping, and especially the linearity of 
the response - that is, the degree to which the response -function amplitudes are propor- 
tional to the magnitude of the input coefficient k. The response functions or shock 
motions due to the application of velocity, mass flow, pressure, or choking-area half-sine- 
wave transients are presented in figures 6, 7, 8, and 9. Examination of the response func- 
tions of the solid-wall inlet (fig. 6(a)) shows that those which are under control (that is, 
the shock position returned to its initial value) are reasonably similar in shape. The 
period is observed to increase from 8 units for k = -0.050 to 14 units for k = -0.100. 
This increase in period appears not only in this case but also for the perforated-wall 
cases. 


It may also be observed that the position to which the shock finally settles is some- 
what different from the initial position of the shock. This change of shock position for the 
same end conditions is believed to be due to a computing phenomenon as a similar change 
occurs regardless of which end condition is held steady (fig. 10(a)). The phenomenon is 
thought to occur because as the computation of the steady-state end condition progresses, 
the locations and number of the points on the characteristic line change. Such changes 
can be expected to cause minor changes in the velocity distribution, and hence a slight 
change in the final steady- state velocity distribution and in the final shock position. 

The effects of opening the wall 20 percent in the region between the shock wave and 
station 4.4 are shown in figure 6(b). The wall perforations increase the value of the maxi- 
mum amplitude of the transient which can be returned to steady state from 0.10 to about 
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0.30, and thus permit a much greater range of controlled operation of the inlet. Also the 
damping of the amplitude of the shock motion has been reduced, as evidenced by the fact 
that the shock overshoots in the downstream direction on its return to its final steady- 
state position. The shock is also returned from a position somewhat upstream of the 
throat, that is, from about station 4.0 to a point beyond station 4.24. The comments con- 
cerning linearity and change of period for the solid-wall case apply generally for the 
20-percent open walls. 

For the transient which was uncontrolled (an amplitude of 0.31), it was observed 
that a secondary shock was formed downstream of the throat with no characteristic points 
existing between the two shocks. This large range between computation values could be 
the reason that the shock does not return to its original position, since all the intervening 
effects of the inlet dimensions are missing from the computation. 

The effects of adding more perforation to the wall, making it 40-percent open to 
station 4.4, are shown in figure 6(c). Significant noise is observed in the controlled state. 
Also, the maximum amplitude of input transient that can be controlled has been increased 
from about 0.30 to 0.50, and the shock returned from a maximum upstream location of 
about 3.7 units. The first uncontrolled run, an amplitude of 0.60, shows a rather violent 
oscillation and little damping compared with the controlled cases. 

The results of the computation using a velocity transient at the engine face show that 
bleeding the subsonic flow in the region between the shock and station 4.4 will allow an 
appreciable increase in the magnitude of the disturbance that can be controlled. The 
appearance of the curves also indicates that linear transfer functions can be used to 
represent the relation between a velocity engine-face transient and the shock location as 
long as the shock location is controlled and the change in period due to larger intensity 
disturbances is taken into account. 

Similar phenomena are observed if either the mass flow rate or the pressure rather 
than the engine-face velocity is allowed to take the transient variation. The effects of 
these transients are shown in figure 7 for the mass-flow variations and in figure 8 for 
pressure variations at the engine face. As with the velocity transient, it is observed that 
increasing the wall porosity increases the amplitude of the transient that is controllable, 
from some value less than 0.005 for the solid wall to 0.200 for 20-percent-open and 
0.250 for 40-percent-open walls in the case of a mass-flow transient, and correspondingly 
from 0.050 to 0.175 for the pressure transient. The observation that the shock motion 
apparently is not controlled, no matter how small the transient, with zero porosity is 
probably due to the previously discussed fact that the initially calculated steady state is 
somewhat different from the finally calculated state and is therefore a calculation phenom- 
enon rather than a physical one. Since the mass-flow value to which the flow must return 
after the transient has passed is set by the initial steady-state calculation, its value is 
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different from a true steady- state final value. This difference in mass flow is compen- 
sated for only by motion of the shock. Hence the shock does not stand still, as would be 
expected, but moves toward the nose to compensate for the changed end condition. 

These unexpected effects of mass-flow transients on the shock location were recon- 
sidered by using another method of controlling mass flow at the engine face. This method 
consisted of restricting the area at the engine face to such a value that the flow through 
the engine face would be choked. The time transient could then be applied to the choking 
area, thereby changing the mass flow in a fashion corresponding to experimental methods 
of changing mass flow. Several steady-state runs were required before a value of the 
choking area A c h (see fig. 1(c)) could be selected which would result in an initial shock 
position reasonably close to 4.4 units. Several of these runs are shown in figure 9(b). It 
is seen that if A c j 1 is chosen to be equal to 0.284 the shock position settles to a value of 
4.467 after a time of about 14. In order to alleviate the effects of the flow adjustment to 
final steady- state values, the transients were not started until a time of about 12 at the 
engine face or about 14 at the shock. The results of applying a half-sine-wave transient 
to vary the choked area at the engine face in a solid-wall tunnel are shown in figure 10(a). 
It is observed that the transient response is very similar to the response observed for the 
velocity change (fig. 6(a)). Also, the shock location returns to the value it had at the start 
of the transient. Thus, none of the irregular shock motions induced by the mass-flow 
transient are observed. 

The effects of applying the half- sine-wave transient to the inlet with distributed 
perforated wall are shown in figure 10(b). It is seen that the magnitude of the transient 
for controlled shock motion is extended from about 0.03 to about 0.3. Also, the controlled 
shock-motion responses are very similar to those observed for the velocity or pressure 
transients. Again the shock position returns to its original value. However, for k = 0.4 
the shock does not go out of control but exhibits a new oscillation and a curve of very dif- 
ferent shape, indicating the possible introduction of a limit cycle if the shock is held too 
long in the region ahead of the throat. 


Shock Motion Due to Engine Stall 

The disturbances to the inlet flow due to engine stalling were simulated by allowing 
the flow velocity to vary as 


u = u n 


i + £(i 

2 


cos o>r) 


( 8 ) 


This function reduces periodically the final velocity of the flow in a fashion very similar 
to that produced by stalling the flow at the engine compressor. Similar disturbances 
could also be computed for pressure, mass flow, or choked end-area variations. The 
results of several calculations using equation (8) are presented in figures 11 to 15. 
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The results of varying the length of the stall from 1 cycle to 6 cycles, with a period 
of 5 units per cycle or about 60 cycles per second, in an inlet with an entrance height of 
1 m are seen in figure 11. For these data, the inlet was 40-percent open from station 0 
to station 4.2 and k = -0.4. Thus the velocity was dropped to 0.6 its original value at 
the minimum value of each cycle. Several phenomena may be observed that the shock 
goes no farther forward than station 3.85. Thus, the shock motion is well controlled. 
Second, the shock motion has a number of frequencies other than the initially imposed 
transient frequency, as may be seen by the highly irregular motion of the shock which is 
indicative of strongly nonlinear phenomena. Third, it is seen that the shock returns to 
its final steady- state end position within a period or two after the stall is stopped, thus 
indicating that the subsonic porous wall is adequately controlling the stalled flow. 

For comparison, figure 12(a) shows the effect of the stall transient on the shock if 
the wall is solid. It is seen that the shock is not contained when the amplitude of the dis- 
turbance becomes greater than about k = -0.10, whereas the perforated walls control the 
shock motion for values of k equal to -0.4 in figure 11 and even higher, -0.5 and -0.6, 
in figure 12(b). 

A test of the linearity (with k) of these flows is seen in figure 12(b), which shows 
the results of using k values of -0.5 and -0.6 for 8 cycles. The shapes are generally 
similar with regard to the location of maximum and minimum amplitudes; however, one 
curve is not a multiple of the other. 

An imposed limit cycle, as observed in figures 9 to 12, is seen more clearly in fig- 
ure 13, which presents a long-period stall (50 time units). It is seen here that a cycle 12 
to 18 units long practically controls the nature of the flow. In fact, the effects of the 
imposed 50-unit cycle can scarcely be seen. It is observed, though, that with the longer- 
period transient cycle the flow returns to normal or steady- state values just as quickly 
as with the shorter-period transient. 

Several stall- simulation runs were made at periods of 12 to 18 units and intensities 
of -0.20 per cycle. The results of these computations are shown in figure 14. The shock 
response is about the same for all the periods tested. Although a slight maximum is 
observed at the period of 16 units per cycle, the increase is so slight that no strong reso- 
nant phenomena can be claimed for this inlet. Also, as in previous cases the oscillatory 
shock motion damps out quickly once the forcing function is stopped. 

The results of several stall-simulation variations of the choked area at the engine 
face are shown in figure 15. This figure also shows the formation of a strong oscillation, 
with weaker shocks that appear to be due to intersection of Q- lines. These weaker shocks 
appear in most of the runs that have the larger amplitude disturbances. In most cases 
they were not plotted so as to reduce confusion in the figures and emphasize the normal- 
shock motion. 
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Control of Throttle -Chop Phenomena 

A change of the throttle setting of an engine causes a transient at the engine face 
such that the final values of the flow are changed from one steady-state value to another. 

It is possible to simulate such a transient by causing one of the flow properties u e , m e , 
or p e at the engine face to experience a ramp-type transient which may be expressed 
as a linear change of ue ? m e , or p e with respect to time from the original steady- 
state flow value to a new steady- state value that is (1 + k) times the original value. The 
transient used consisted of a linear portion extending over a time of 40 units, or about 
one-eighth of a second in an inlet of 1- meter height, after which the flow retained a new 
constant value. 

A variety of cases were investigated. Several velocity ramp transients were run 
with a solid-wall nozzle in order to establish a basis for comparison. (See fig. 16.) The 
perforation distribution presented in figure 5 was used for a number of runs with varia- 
tions in either u e or the choked area at the engine face. The results of these runs are 
presented in figures 17 to 19. Figure 16 shows that with solid walls the shock position is 
bounded for only a very limited change in the value of the engine-face velocity. As would 
be expected, once the shock is pushed ahead of the throat it moves immediately out of the 
inlet. This phenomenon occurs for a reduction in velocity of about 3 percent. The addi- 
tion of the distributed porous wall shown in figure 5 causes remarkable changes in the 
shock motion. In figure 17 it is seen that reductions of the engine-face velocity by as 
much as 60 percent (k = -0.6) do not cause the shock to leave the interior of the inlet. 

Even in an extended run, to 200 units rather than 80 units, the shock is contained (fig. 18) 
to some value greater than r = 200, or about 0.6 sec in an inlet 1 meter high. The appar- 
ent escape of the shock may be due to an accumulation of errors because of the length of 
the computer run required. However, in either case, the perforated wall used in conjunc- 
tion with other inlet control systems does allow sufficient time for these systems, such 
as area distribution change or a bypass door opening, to return the engine-face operating 
conditions to their original steady- state value. 

The question of the existence of this oscillating phenomenon in an actual porous- 
wall inlet cannot be decided, as the possibility exists that an unstable computing condition 
is inherent in the computing program. However, in any investigation of this type of inlet 
it would be well to include instrumentation capable of checking for such oscillations. It 
is believed probable that whenever the shock is forced ahead of the throat such oscilla- 
tions can be expected. For example, in figure 9(a), the oscillations are seen to develop 
even though the end conditions are constant if the shock is started ahead of the throat, 
whereas no oscillation occurs if the shock is started downstream of the throat. In fig- 
ure 9(b), with A c h = 0.324 and A c h = 0.284 the shock position shows no sign of oscil- 
lation. It is believed that if computing instabilities existed they would also be evident in 
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these cases where the shock remains downstream of the throat. The calculated oscilla- 
tion is therefore believed to give a valid indication of the behavior of an actual perforated 
supersonic inlet under a throttle- chop velocity disturbance. 

Similar phenomena are shown in figure 19, which presents the results of several 
computations using a ramp variation of the choked area at the engine face. It is seen 
that the distributed perforated wall controlled the flow for a change in of as much 

as 0.33. It is also seen that the choked-flow end condition shows the limit cycle strongly 
when the shock location moves ahead of the throat. 

CONCLUDING REMARKS 

This paper presents an analysis of the ability of a perforated wall in an internal- 
compression supersonic inlet to control the motion of the shock due to a variety of tran- 
sients in the flow values at the engine face. The following remarks are subject to the 
possible limitations imposed by computing with a finite-increment one-dimensional char- 
acteristic system. 

1. The use of internal-compression inlets with walls selectively perforated for 
subsonic flow makes possible containment of the shock within the inlet for much larger 
engine-face transients than with solid walls. 

2. Selectively perforated walls contain the shock ahead of the throat for reduced 
flow properties at the engine face; however, a limit cycle of roughly constant magnitude 
is set up by reflection of the waves between the shock surface and the engine face. 

3. Some approximation to linearity of responses to impulse-type transients was 
observed. However, the development of the limit cycle eliminates any appearance of 
linearity in the responses to longer (nonimpulse) transients when the average shock 
position is ahead of the inlet throat. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., May 2, 1969, 

720-03-00-02-23. 
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(a) Design of the perforated wall. 
Figure L- Inlet design. 
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Subsonic pressure opens flapper valves 
Springs close valves for supersonic flow 

(b) Structure of perforated wall considered for control of shock position. 
Figure 1.- Continued. 
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(a) Solid wall. 

Figure 6.- Shock motion due to sinusoidal impulse in engine-face velocity. 
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